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ABSTRACT 

Scour About A Single, Cylindrical 
Pile Due to Random Waves and Current. (December 1984) 
Robert W. Eadie IV, B.S., United States Naval Academy 
Chairman of Advisory Committee: Dr. J. B. Herbich 

An empirical approach was utilized in this investiga- 
tion to study the scouring effects that various currents 
and random waves have on a single cylindrical pile. Re- 
sults from laboratory experiments are presented in graphi- 
cal form depicting relationships between pertinent dimen- 
sionless parameters which were developed through a dimen- 
sional analysis. These results were then compared to re- 
sults compiled in previous studies that investigated scour 
due to currents and monochromatic waves in order to deter- 
mine if the irregularity of the wave train is significant 
in the study of scour. 
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PART I 

INTRODUCTION 

The construction of offshore structures is an ex- 
tremely expensive process. A portion of these construc- 
tion costs are due to large safety factors used during 
the design process. These safety factors which may lead 
to overdesign and, therefore, unnecessary expense are the 
result of the engineer's inability to accurately predict 
some of the forces on the structure, as well as some of 
the natural phenomena that will occur once the structure 
is in place. One of these natural phenomena which is be- 
coming an ever increasing problem is scour caused by waves 
and currents around offshore structures. Design factors 
for offshore structures depend on the location of the ocean 
bottom. The scour of a non-cohesive ocean floor can lower 
the bottom line significantly, and by doing so can dras- 
tically change these design factors. 

The purpose of this investigation was to study what 
effect various currents and random waves have on a singu- 
lar, cylindrical pile. These results were then compared 
to results compiled in previous studies that investigated 
scour due to currents and monochromatic waves to determine 
if the irregularity of the waves is significant. An em- 
pirical approach was used in this study. 
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Experiments were conducted in a two-dimensional wave 
tank. The pile used in this study had a diameter of 1.5 
inches. Two water depths, two currents, one sediment size 
and four random wave spectra were utilized. Using data 
obtained from the experiments, an attempt was made to de- 
scribe scour in terms of relevant dimensionless and dimen- 
sional parameters. By doing so, it is intended that the 
reader will obtain a better understanding of the scouring 
process . 
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PART II 

REVIEW OF LITERATURE 
Incipient Motion 

The scour phenomenon begins with the incipient motion 
of sediment particles. Eagleson and Dean (7) define in- 
cipient sediment motion to be "an instantaneous condition 
reached when the resultant of all the active forces on 
the particle intersects the line connecting the bed parti- 
cle contact points." There are three types of active 
forces acting on a sediment particle. Two of these forces, 
the drag and lift forces-, are due to the fluid motion. 

These forces are in turn balanced by the sediment parti- 
cle mass. These forces are depicted in Figure 1 (1). 

A thorough understanding of these forces is required 
in order to comprehend fully the processes of incipient 
motion, bed-load transport and bed form development; all 
which depend on the collective response of initially sta- 
tionary bed particles to fluid forces. The following dis- 
cussion of incipient motion is for a unidirectional flow. 

The variation of drag forces with flow condition has 
been well documented, as has the lift force in hydrodyna- 
mically-rough flow conditions. However, investigations 
of lift forces in hydrodynamically smooth and transitional 
flow have yielded conflicting results. 
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Fig* 1 . 



Hydrodynamic force on an individual sediment 
particle (1 ) 



Prior to 1967, most investigations measuring fluid 
lift forces were conducted in the hydrodynamic rough flow 
region characterized by: 

Ro. = Vj.D > 70 

A /V 

V 

in which IL. represents the grain Reynolds number; Vj f = 
shear velocity; D = grain diameter; and v = kinematic vis- 
cosity of the fluid. All measurements under these hydro- 
dynamic conditions indicated positive lift. 

Coleman (4) was the first to study lift forces in 
the hydrodynamically smooth region (R* < 5) and transition- 
al flows (5 < R* < 70) . His work reported measured drag 
forces and inferred lift forces from observing the inci- 
pient motion of spheres having various densities. From 
the results, it was suggested that for a sphere resting on 
top of a closely packed bed of spheres of equal diameter, 
the lift force decreased through zero and became negative 
(bedward-directed) as Rj, decreased through a value of ap- 
proximately 15. 

Davies (6) prepared an experiment whose sole purpose 
was to define the flow conditions under which negative 
lift forces occur. He found that significant negative 
lift can be expected to occur only if the following con- 
ditions are met: 1) The value of R.* < 5 and 2) The test 

sphere is situated such that a significant volume rate of 
flow can occur beneath the sphere. In a situation in 
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which considerable fluctuation of bed shear stress occurs, 
the onset of significant negative lift with a decrease in 
R*> will occur gradually and unsteadily as the effect of 
turbulence in the area of the bed particle becomes less 
significant, until hydrodynamically smooth flow is es- 
tablished at Rj. < 5 and the lift force becomes predominant- 
ly negative. 

The existence of varying types of lift forces is im- 
portant with regard to bedform development. If the par- 
ticle is removed from its position of rest at R v . < 5, the 
negative lift forces acting on the particle will tend to 
prevent the particle from leaving the immediate vicinity 
of the bed, and, therefore, it will move along the bed as 
a bed load. On the other hand, if the particle is removed 
at R* » 5, it will experience a positive lift force that 

will tend to impart an upward acceleration and thereby 
lift the particle into the main turbulent flow, where it 
may or may not become suspended. 

There is considerable evidence to suggest that a rip- 
pled bedform exists only at Rj. < 5, and that a dune bedform 
exists on its own only at R* > 24. It is therefore sug- 
gested that: 1) The flow conditions that are necessary 

for existence of the negative lift forces on a sediment 
particle correspond closely to those required for a plane 
or rippled bedform. 2) Positive lift force conditions 
correspond to a dune bedform. 
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The condition of incipient motion for a cohesionless 
soil particle is described in terms of the forces acting 
on the particle by: 

F^ 



tan 0 = p— 
n 



( 1 ) 



where F and F^ are the forces parallel and normal to the 
angle of repose 0. In our case, F^_ and F^ are the resultant 
forces of the three active forces mentioned earlier, the 
hydrodynamic drag force F^, the lift force F^ , and the 
submerged weight W (gravity force) . In accordance with 
equation (1) , the condition of incipient motion under the 
influence of these three forces becomes: 

Ws ina + F , 



tan 0 = n— w- 

Wco sa - F, 



( 2 ) 



where the angle a is the inclination of the bed from the 
horizontal at which the incipient motion occurs. The above 



situation is illustrated in Figure 2. 

2 



F D = C D K l d 



2 pU b 



(3) 



and 



F L = C L K 2 d ' 



pU. 



( 4 ) 



where 

U, = bottom fluid velocity 
b 

= drag coefficient 

C T = lift coefficient 
.Li 



8 







Fig. 2. - Force diagram on particles in a cohesionless 
loose bed (9) 
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d = particle diameter 

K 1’ K 2 = s ^ a P e fsctor 
p = fluid density 

The submerged unit weight of the particle is expressed 
as : 

W = k 3 (p s - P)gd 3 (5) 

where 

= shape factor 
p s = particle density 

By introducing equations (3) , (4) and (5) into equation 

(2), the result becomes, 

2 

(U^ ) cr 2 K 2 (tan 0 cos a - sin a) 

( p s / p -l)gd = c d k i + C l K 2 tan 9 (6) 

where (U^)cr is the critical bottom velocity at which ac- 
cording to equation (1) incipient motion occurs. 

(U b )cr =V T o/ p (7) 

where 

t = critical shear stress 
o 

p = fluid density 

The value of the right hand side of equation (6) de- 
pends on A) The particles: their size, shape, size dis- 
tribution, uniformity, . . . ; B) The flow dynamics, which 

determine the values of and C^; C) the channel slope 
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which in many cases amounts to cos a ~ 1 ; and D) The angle 
of repose which depends on the particle properties (9) . 

Once the critical flow velocity is reached, the cri- 
tical shear stress will also be obtained as per equation 
(7) . Most data on critical shear stress for non-cohesive 
sediments have been observed in flume experiments. These 
experiments have shown that the motion of the sediment 
particles at the bed for uni-directional flow is highly 
unsteady and nonuniformly distributed over the entire bed 
area. Near critical conditions, the particle motion in 
any small area of the bed occurs in gusts whose incidence 
increases as the mean shear stress increases. When ob- 
serving a large section of a sediment bed as the shear 
stress approaches the critical value, the occurrance of 
these gusts of particle motion seem to be random in both 
time and space. 

Two of the best known and widely accepted investiga- 
tions into the concept of critical shear stress are those 
of Shields and White (17) . Shields developed the relation: 




Uj,d 



■) 



( 8 ) 



v 



where 



t = critical shear stress 



o 



Y = specific weight of sediment 
s 

Y = specific weight of fluid 
d = particle diameter 



11 



U. v = shear velocity 
v = kinematic viscosity 

Shields' results, in which he determined the function 
U*d 

f (— — ) is shown in Figure 3. 

As discussed earlier under Davies experiment, under 
flow conditions in which R v , < 5, the particle will have 
negative lift and will therefore be less likely to move. 
This corresponds to the Shields curve that indicates that 
when Rjj. decreases below about 10, particles increasingly 
resist the dimensionless shearing stress of a flow. 

White (17) determined the critical shear stress for 
several sediments under conditions of developed laminar 
flow, laminar boundary layer flow and turbulent boundary 
layer flow. White concluded that the true critical shear 
stress necessary to move a sediment particle has a fixed 
value and that this value is given by experiments dealing 
with laminar flow. White's (24) equation for t c for sedi- 
ments in a horizontal bed is 

t = 0.18 (v - y)d tan 0 (9) 

c ' s 1 s 

where 

2 

t = critical shear stress (dynes/cm ) 
c 

0 = the angle of repose of the sediment submerged 
in the fluid 

and in which the constant (0.18) is obtained from experi- 
ments with laminar flow. 



Dimensionless shear stress, T 



12 




4 6 8 10 20 40 

Boundary Reynolds Number, R , = 



Fig. 3 » - Shields diagram with White data added ( 24 ) 
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Yang (28) expressed some reservations on using Shields' 
diagram as the criterion for incipient motion. Using basic 
concepts in fluid mechanics, the boundary layer theory 
and laboratory and field data from various investigators, 
Yang plotted the relationship between dimensionless criti- 
cal velocity and shear velocity Reynolds number at inci- 
pient motion (Figure 4) . From this curve, the following 
equations of incipient motion were obtained: 



V 



CR 

W 



2.5 

raid 



U . d 



+ 0.66; for 0 < 



< 70 . 



log (-4L-) - 0.06 



^CR 

-T-2.05; for 70 < — 



( 10 ) 

(ID 



where 

Vp V — average critical velocity at incipient motion 

U Xa. 

W = terminal fall velocity, which is reached when 
the particle drag force and the submerged unit 
weight of the particle are balanced, 
d = particle diameter 
v = kinematic viscosity 

Uj. = shear velocity (Uj. = '^ T / , P = fluid density, 

t = average shear stress) 

Uj.d 

— — — = shear velocity Reynolds number 

v 

Many experiments have shown that for any given bed 
material, there is no single value of bed shear stress or 
velocity, below which not a single particle will move and 



oixniMim cmticm. vttocm 
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Fig* 4 



Relationship between dimensionless critical 
velocity and shear velocity Reynolds number 
at incipient motion ( 28 ) 



L J 



above which all the particles having the same size will 
move. The forces acting on the particle are not constant, 
but instead fluctuate about some average value, a phenomenon 
closely associated with turbulence. Also, the forces which 
move the particle are dependent on the particle position. 
Therefore, as mentioned earlier, the movement of bed ma- 
terial particles is very unsteady, and non-uniformly dis- 
tributed over the entire bed (17) . Einstein (8) has sug- 
gested that the process of bed material movement is sta- 
tistical in nature. Because of this, it appears that no 
truly critical condition for particle movement exists. 

This is further exemplified in Figures 5 and 6 which show 
the variance in values for the critical shear stress and 
critical water velocity for quartz sediment. 

Scour 



Fluid Mechanics of Scour 

Local scour is defined as "the abrupt decrease in 
bed elevation near a pier due to erosion of bed material 
by the local flow structure induced by the pier." (22) 

One type of classification of the type of scour is the 
amount of sediment transport into and out of the scour hole. 

Laursen (12) suggested that the scour process could 
be represented by the relationship: 

dQ s 

at - " q s i “ q s2 



. ( 33 ) 



Critical shear stress, T c , In pounds per square feet 




Fig. 5. - Critical shear stress for quartz sediment in 
water as a function of grain size (24) 




Fig. 6. 



Critical water velocities for quartz sediment as 
a function of mean grain size (24) 



Velocity, in centimeters per second 




in which dQ^/dt — the rate of local scour in volume per 
unit time; q gl = the capacity of the flow to transport 
sediment out of a scour hole in volume per unit time; and 
q s 2 = the rate at which the undisturbed flow supplies sedi- 
ment to the scour hole. 

From this expression, three cases of scour may occur. 

They include: 1) no scour, in which case q n = q 0 = 0; 

si s Z 

2) clear water scour, where q > 0 and 0 ~ q 9 « q ; and 

Si * s z si 

3) scour with continuous sediment motion, where q , > q 0 > 0. 

n sl n s2 

Abrupt changes in bed elevation near a pile occur only in 
cases (2) and (3) . These two cases of scour are depicted 
in Figure 7. 

The main flow feature responsible for scour whenever a 
slowly-varying boundary layer current is present, results 
from secondary flow. The passage of a steady flow boundary 
layer along the surface past the cylinder causes a stagna- 
tion pressure gradient. Because the stagnation gradient 
decreases 'toward the sediment bed, a secondary flow is 
generated in the direction of the bed. This secondary flow 
has two effects on the shear or boundary layer flow. The 
principle effect is the formation of a horseshoe vortex 
that wraps itself about the cylinder and trails off down- 
stream. The horseshoe vortex shown in Figure 8 is depicted 
as a spiral curve extending around the cylinder base (15) . 

Considerable laboratory data are available from scale 
model studies; however, most of these studies investigate 



Scour d*pth, dg 
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M (b) 



Fig. 7* - Time variation of scour depth (22) 
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Fig. 3 . - Secondary flow about a cylinder ( 15 ) 



the horseshoe vortex in a river environment. Even though 
these results provide some insight into the characteris- 
tics of flow and the basic mechanism of generation of the 
horseshoe vortex, it is impossible to apply these river 
model results directly to an ocean environment. A major 
deficiency being the lack of appropriate dynamic and kine- 
matic scaling laws. Shen et al . (21) have shown that the 

size and the strength of the horseshoe vortex can be re- 
lated to a Reynolds number using the characteristic length 
of the cylinder in the flow field. This relationship is 
as follows: 

a 2aU 

Vortex strength = = f( — — ) (34) 

where 

u = rotational velocity of vortex core 

A = area of vortex core 

a = half width of pile 

U = undisturbed flow used in the potential flow sense 

oo 

v = kinematic viscosity of fluid 

Also, Shen et al. (21) suggested that for cylinders 
with a height greater than the thickness of the bottom 
boundary layer, the area affected by the horseshoe vortex 
extends approximately one cylinder diameter upstream of 
the obstruction. 

Hjorth (10) indicated that the fluid shear stress 
striking the bottom below the horseshoe vortex system are 



increased by a factor of eleven when compared with the 
shear stress values associated with the undisturbed free- 
stream. This increase in bottom shear stress causes the 
bottom sediments to scour. With scour occurring, the bot- 
tom elevation is altered, which in turn brings about a 
change in the flow in the horseshoe vortex. 

Another effect of the secondary flow about a cylindri- 
cal obstruction is the downward deflection of the external 
flow streamlines about the cylinder. This downward deflec- 
tion of streamlines is shown on Figure 8, the secondary 
flow is depicted at the stagnation line on this same fig- 
ure. The convective effect of the streamlines has little 
effect on the horseshoe vortex, except for a slight strength- 
ening of the rotational vortex. 

A circular cylinder blocks the approach flow, result- 
ing in convective acceleration. The acceleration causes 
an increase in the bottom shear stress by a factor of up 
to four over the bottom shear stress of an unobstructed 
flow. 

Flow patterns in the cylinder wake also influence 
scour. There is a net downward flow in the wake of a 
cylinder resting on a plane surface. This is due to the 
fact that the wake region is a low pressure region. The 
wake in the high velocity region is at a higher pressure 
than the wake in the low velocity region near the sedi- 
ment bed. This difference in wake pressure causes a flow 



to form in the wake acting downwards toward the bed. If 
the effects of the horseshoe vortex are included, this 
downward flow is enhanced ( 15 ) . 

In an attempt to get a better understanding of the 
mechanisms of scour, the subject will be dealt with in 
three parts. First, the mechanics of scour in a steady 
current alone will be discussed, followed by the mechanics 
of scour due to oscillatory flow, and, lastly, the mechanic 
of scour caused by the interaction of waves and currents. 
Mechanics of Scour in Steady Currents 

Repeated observations indicate that the scale and 
shape of the scour hole closely resemble the scale and 
shape of the horseshoe vortex. The scour is caused by 
the large bottom fluid shear stress gradients occurring 
beneath the horseshoe vortex. The initial vortex, before 
the onset of scour, exerts bottom fluid shear stresses 
that are approximately one order of magnitude greater than 
the shear stresses in the surrounding undisturbed flow. 
Therefore, it is possible for scour to occur even though 
sediment is not eroded and transported into surrounding 
regions. If general bedload transport is occurring, then 
the elevated bottom shear stresses beneath the horseshoe 
vortex produce more rapid sediment transport and bed ero- 
sion . 

As the erosion of sediment continues, the intensity 
of the horseshoe vortex decreases as it spreads out into 



the scour cavity. As the vortex intensity decreases, the 
fluid shear stress acting on the sediments in the bottom 
of the scour cavity are also decreased. Once equilibrium 
is reached, the fluid shear stress acting upstream at the 
bottom of the scour cavity is equal to the fluid shear 
stress acting downstream on the undisturbed portion of 
the bed. The fluid shear stress acting on the upstream 
face of the scour hole is equal to the fluid shear stress 
acting on the undisturbed area of the sediment bed, plus 
an increment approximately equal to the fluid shear stress 
required to entrain the sediment. Once this distribution 
of shear stresses has been achieved, the scour hole has 
reached its equilibrium size and shape (Figure 9) and the 
rate of transport of sediment from the scour hole equals 
the rate of supply (15) . 

Mechanics of Scour for Oscillating Flow 

Observations of scour development from oscillating 
flow are given by Ranee (19) and Palmer (18) . The exact 
sequence of scour hole development depends on the diameter 
of the cylinder, the intensity of the oscillating flow, 
and the sediment transport on the surrounding bed. 

Ranee states that for a slender body (i.e., the dia- 
meter of the body is small when compared with the wave- 
length) under the action of waves, the shape of the scour 
hole is dependent on the amplitude of the wave motion at 
the bed. When the amplitude is very small when compared 
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INITIAL STAGE 




Fig. 9. - Steady current occur cavity (15) 



Co pile diameter , two discrete scour holes are formed on 
the forward part of the cylinder at 45° to each side of 
the centerline. For this case, eddy shedding has no ef- 
fect on sediment transport. When the bed semi-orbit length 
is larger and becomes approximately one fifth of the pile 
diameter, there is a tendency for flow separation to occur. 
Accretion and erosion are localized and four scour holes 
versus two are formed. These scour holes have a maximum 
depth of about 10% of the pile diameter. At larger ratios 
of semi-orbit length divicLed by pile diameter, of about 
unity and above, eddy formation is pronounced and a scour 
hole develops completely and uniformly about the pile with 
a maximum depth of 40% of the pile diameter. 

Palmer noted that the scour hole initiates at the 
sides of the cylinder and, in many cases, some sediment 
deposition was found to occur against the upstream and down- 
stream sides of the cylinder. As the scour hole continues 
to develop, these deposits are removed and a majority of 
the sediment deposition shifts to a rim outside the scour 
hole. Once fully developed, the scour hole was found to 
be a radially symmetric frustrum of an inverted cone hav- 
ing side slope at the angle of repose of the sediment. 
Palmer's data indicated that the depth of the equilibrium 
scour hole to the pile diameter ranged from between 0.5 
and 0.75, a value higher than that previously stated by 



Ranee . 



For oscillating flow past a cylindrical obstruction, 
the horseshoe vortex is insignificant, and, therefore, 
scour is initiated from the bottom fluid shear stress 
caused by the acceleration of primary flow past the ob- 
struction. This results from a bottom fluid shear stress 
gradient that is positive from upstream of the cylinder 
to the point of maximum velocity, and negative from that 
point to the downstream limit of undisturbed flow (Figure 
10) . If the maximum bottom fluid shear stresses are great- 
er than the sediment entrainment threshold, erosion will 
occur from upstream of the cylinder up to the point of 
maximum velocity, and deposition will occur downstream of 
this point. 

In order for an equilibrium scour hole to develop 
around the entire circumference of the cylinder, an addi- 
tional flow condition must develop. The initial scour hole 
produces bottom slopes close to the pile which must even- 
tually cause flow separation in the thin orbital boundary 
layer. Reverse eddies in the separated flow erode the 
initial side slope of the scour hole until it obtains a 
value of the angle of repose of the sediment. With the 
steepening of the sideslopes, the zone of separated bottom 
boundary layer flow spreads out around the pile eventually 
causing the scour hole to obtain the equilibrium shape 
shown in Figure 11. Therefore, the equilibrium scour hole 
for oscillatory flow is related to boundary layer separation 
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Fig. 10. - Sediment transport due to oscillatory flow (15) 
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INITIAL STAGE 





EQUILIBRIUM STAGE 



Fig. 11. - Oscillatory flow scour cavity ( 15 ) 



and not dependent on the nature of secondary flow as was 
the case for scour due to steady currents (15) . 

Mechanics of Scour Due to Waves and Currents 

In this study, the scour that developed through the 
combined action of waves and currents was investigated. 

This is the type of scour that occurs in almost all ocean 
environments. Ninomiya et al. (16) performed scale model 
experiments using the combined action of waves and currents. 
They found the shape of the scour hole to be similar in 
shape to that developed by a steady current without waves, 
and that the rate of scour hole development was consider- 
ably faster in the case of combined wave and current action 
when compared to steady current. 

As mentioned earlier, if the steady flow velocity is 
stronger than the oscillating flow velocities, a horseshoe 
vortex will form about the base of the pile. The intensity 
of this vortex varies according to differences in the 
square of the -combined steady and oscillating flow velo- 
cities. The size and shape of the equilibrium scour hole 
should not be as large as that caused by a steady current 
whose mean velocity is equal to the combined velocity of 
the wave orbital and mean current condition. This is due 
to the fact that velocity profile of the combined veloci- 
ties will be flatter than that achieved in a simple 
steady boundary layer flow. On the other hand, the equi- 
librium scour hole depth should be greater than that caused 
by a steady boundary layer having no wave action. 



If a strong current and a weak wave orbital flow 
exists, the development of the scour hole will depend 
primarily on a horseshoe vortex with a pulsing intensity. 
The scour hole development will be faster when compared 
to a steady current alone due to the additional bottom 
fluid shear stress caused by the wave orbital boundary 
layer and because the vortex formed from the steady por- 
tion of the flow can more easily transport the wave- 
agitated sediment. 

In the case of a strong orbital flow and a weak cur- 
rent, the scour hole shape will be similar to that caused 
by a steady current. Once again, the scour will develop 
more rapidly than it would under the influence of a steady 
current alone due to the additional bottom fluid stresses 
caused by the waves orbital motion. In this case, the 
equilibrium scour depth can be less than that associated 
with a steady current that is equal to the sum of the weak 
current and orbital velocity. 

The situation becomes difficult to evaluate when the 
magnitude of the oscillatory motion and the steady current 
are nearly equal. The scour hole size and shape will nor- 
mally depend on the horseshoe vortex associated with a 
steady current. It is not possible to generalize on the 
scale and intensity of the pulsing vortex when compared 
to the steady current situation. 






The size and depth of scour produced by the combination 
of waves and currents acting on a pile that is large in 
comparison to the oscillation amplitude will be drastical- 
ly smaller than that predicted for a small pile, except 
when the steady current is greater than the maximum wave 
orbital velocity. This is caused by the fact that the 
horseshoe vortex does not have time to wrap completely 
around the pile before the flow reverses due to the wave 
orbital velocity (15) . 

Previous Studies on Scour 

Very little experimental work has been done on scour 
caused by the combined effect of oscillatory wave motion 
and current; however, there have been numerous studies 
dealing with scour in open channel flow. Because the 
scouring forces are similar for oscillatory flow as for 
steady state open channel flow, that is, hydrodynamic in 
nature, the results found in experiments in open channel 
flow can be applied with reservation to the combined oscil- 
latory and steady current flow. 

Some experimental results are presented by Shen et 
al. (22) as shown in figures 12-14. These experiments 
employed a steady current. Figure 12 indicates that the 
depth of scour is a function of the pier Reynolds number 

OR) • 

TR = UD/v 



where 
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Scour depth as a function of pile Reynolds number (22) 



Equilibrium scour depth, d $a , feet 
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Fi^ 13* - Equilibrium scour depth - Pier Reynolds 
’ number relationship with various pier 
sizes (15) 
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Fig. 14. - Equilibrium scour depth - Pier 
Reynolds number relation with 
different sediment sizes (22) 



U = horizontal free stream velocity 
D = pile diameter 
v = kinematic viscosity 

It was also shown that the scour depth is also dependent 
on the pile diameter (Figure 13) and sediment size (Figure 
14) , and this explains why the Reynolds number alone is 
insufficient to collapse all of the data shown in Figure 
12 onto a single line. 

Carstens (3) performed studies of scour associated 

with various types of obstacles. He concluded that the 

rate of scour induced by the obstacle is a function of 

the sediment number (N ) , sediment geometry, as well as 

s 

the ratio of the scour depth to obstacle size. It should 
be noted that the "sediment number" is a densimetric Froude 
number . 



U = free stream velocity 

S = sediment specific gravity 
s 

g = acceleration of gravity 
d = mean diameter of sediment particle 
Coleman (5) used the principle that in any problem 
of flow around an immersed body, where the dependent vari- 
able is a pressure drop, the most general formulation of 
the problem is in terms of an Euler number and either a 




where 



Reynolds number, a Froude number or both, with the dependent 
variable appearing in the Euler number. He derived a rela- 



tion between the scour Euler number (U/^ 2gd ) and the pier 

s 

Froude number (U/ VgD) . 



U 



0 . 6 (-^£> 



U s 9/10 



(35) 






From the above relation and using a previously derived 
depth of scour relation using the pier Froude number, an 
expression was obtained to predict the depth of scour for 
steady flow. 



d g = depth of scour (ft) 

D = pile diameter 
U = free stream velocity 
g = acceleration of gravity 
Wells and Sorensen (27) investigated scour around a 
single circular pile due to oscillatory wave motion. Fig- 
ure 15 shows that a relationship exists between the rela- 
tive scour and the relative water depth for various sedi- 
ments. In addition, their results indicated that the ra- 
tio of ultimate scour to wave height appeared to be di- 
rectly related to the sediment number, N g , (Figure 16) and 
to the pile Reynolds number, Npp , (Figure 17). In Figures 
15 though 17, the notation is defined as follows: 




(36) 



where 



37 




CM 




1 

LA 



to 

•H 



Relative ultimate significant scour depth as a function of 
relative depth (2?) 
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Fig. 16. - Relative ultimate significant scour depth as a function of the 
sediment number (2?) 







Fig. 17. - Relative ultimate significant scour depth as a function of 
the pile Reynolds number (27) 



S u - ultimate significant scour depth 
d = mean grain diameter 

Sand no . 1 : d = . 62 mm 

Sand no. 2: d = .325 mm 

Sand no. 3: d = .30 mm 

h = still water depth 
g = acceleration of gravity 
T = wave period 
H = wave height 

U = maximum free stream velocity 
S = specific gravity of sand 
D = pile diameter 
v = kinematic viscosity of water 
More recent scour studies include the work of 
Raudkivi and Ettema (20) , Jain and Fischer (11) , Armbrust 
(2) , and Wang (26) . The work of Armbrust and Wang dealt 
with the scour associated with the combined effect of waves 
and currents while the remaining studies dealt only with 
a steady current. 

Raudkivi and Ettema (20) performed laboratory experi- 
ments on clear water scour of cohesionless bed sediment 
at cylindrical piles. They found that: 

1. As the geometric standard deviation of the par- 
ticle size distribution increases then the equi- 
librium scour depth decreases. 



2. The equilibrium depth of scour decreases with 
the decreasing relative size of pier diameter 
to mean bed particle size for values less than 
about 20 to 25. 

3. The maximum equilibrium depth of scour is a 
function of the relative depth h/D and the rela- 
tive pier size D/D^q. 

4. The equilibrium depth of scour decreases at a 
greater rate with decreasing flow depth for 
smaller values of the relative flow depth h/D. 

Jain and Fischer (11) performed scour experiments 
for steady currents at high velocities, i.e., for flow 
regimes well beyond the onset of general sediment trans- 
port. From these studies, they concluded that after de- 
creasing slightly upon the start of sediment transport, 
scour depth does increase again with increasing velocitie 
This conflicts with the consensus of previous experiments 
which concluded that once the conditions of general sedi- 
ment motion have been established in a stream, there will 
be no further increase of scour depth with velocity, be- 
cause the rate at which sediment enters the scour region 
is equal to the rate of sediment removal by the scouring 
process . 

Armbrust (2) found that in general, for a given cur- 
rent and relative depth (h/L) , scour increases as the 
wave steepness (H/L) increases (Figure 18) . Figure 18 
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Scour as a function of wave steepness for constant values 
of current and relative depth (2) 



also shows that the ultimate scour depth increases as the 
relative depth decreases. Additionally, in general, for 
a fixed relative depth, Figure 19 reveals an increase of 
scour depth with an increase of relative wave height (H/h) . 
In Figures 18 and 19, the notation is defined as follows: 

= ultimate scour depth 
D = pile diameter 
D 50 = mean sediment diameter 
H = wave height 
L = wave length 
£ = f ree stream velocities 
h = still water depth 

Wang (26) defined ultimate scour (S ) to include 
"general scour" and net "local scour" and, therefore, ul- 
timate scour is equal to the total scour that occurs. Wang 
stated that: 

S u = f (V^ w , h , p , H , L , T , d^Q , p s , D , v , t) 
where 

V c = combined wave and current velocity 
fw 

h = water depth 
p = water density 
H = wave height with current 
L = wave length with current 
X = wave period with current 
D = pile diameter 
d^Q = mean sediment diameter 



Points above dotted line - = .81 fps 

Points below dotted line - U-, = .43 fps 
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Hg« 19 • ~ Scour as a function of relative depth for constant value 
of current and relative wave height (2) 



p = density of sediment 
v = kinematic viscosity of water 
t = time required to produce scour 
Using dimensional analysis, Wang produced the relation: 

S 

= f (N f • H/L • U r • N g • N Rp ) = f (6) 
where 



6 = N-r. ' H/L • U 



N 



N t 



’F “ u r "s RP 
V 3 H 2 L[V + (1/T - V/L) t^-] 2 D 



p s ' p 2 J 

g v d 
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Figure 20 is a plot of the relationship of relative scour 
depth, S^/h, and the parameter 8. 



SAND NO. I 
SAND NO. 2 




PART III 



THEORETICAL ASPECTS 
Wave-Current Interaction 

Hoften and Karaki (23) performed a laboratory study 
to investigate the interaction of waves and currents as 
mechanically generated, monochromatic waves propagated 
on a turbulent open channel flow. 

They found that waves propagating on an open channel 
flow are shown to be modified by a combination of the wave 
interaction with the shear gradient and the existing tur- 
bulence. The effect was to remove energy from the waves 
causing an increase in wave energy dissipation. 

The flow itself is also changed by the waves. The 
induced orbital wave velocities distort the mean velocity 
profile increasing the gradient close to the sediment bed. 
Turbulence energy is formed near the surface by the in- 
teraction of the wave-induced Reynolds stresses and the 
fluctuating velocity gradients. The energy diffuses down- 
ward where it is dissipated on the sediment bed in the 
form of an apparent increase in boundary shear stress, an 
important factor in the scouring process. 

In order to evaluate the combined velocity of a wave 
propagating in the same direction as a uniform current, 
the Airy wave theory will be used. By doing so, if a 
uniform current is assumed, the velocity potential of the 
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uniform current and the velocity potential due to 

wave (4> 2 ) can added directly to obtain the velocity 
potential of the combined wave current condition (<t>) be- 
cause both velocity potentials are linear solutions of the 
Laplace equation. The parameters of wave period (T) , wave 
height (H) and wave length (L) are shown in Figure (21) , 
and are properties of a wave propagating in still water. 
Wave velocities can be calculated using these parameters 
with the water depth for a no current condition. There- 
fore, the velocity relation for combined wave-current in- 
teraction becomes, 



+ <J> £ (12) 

where 

<j>, = U (13) 

T 1 c 



and the derivation of <J> ^ I s outlined in the following pages. 
This procedure closely follows Newman (14) . 

The fluid velocity under wave action is expressed by 
the gradient of the velocity potential <p , and the effects 
of the free surface must be expressed in terms of appropri- 
ate boundary condition on this surface. 

The exact kinematic boundary condition can be derived 
most readily by requiring that the substantial derivative 
of the quantity z - n vanish on the free surface. The 
result of this condition is that on z = n, 

D , v _ 3<£_ Sn _ 9<t> Sn S4> 4>n 

0 “ Dtf ^ z - n ' dz ~ St Sx Sx ay ay 



. . ( 14 ) 
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For the linearized kinematic boundary condition, the second 
order terms can be dropped, which gives the relation: 

3n _ 3<j> 



This approximate boundary condition states that the verti- 
cal velocities at the free surface and the fluid particles 
are equal, disregarding small departures of the surface 
from the horizontal orientation. 

The dynamic boundary condition is obtained from 
Bernoulli's equation. The exact condition to be satis- 
fied on the free surface is : 

‘ k (P ‘ P a ) = ft + 1/2 7 ♦ ’ V * + § z = 0 (16) 

Substituting the free surface elevation n for Z and solv- 
ing for n the expression becomes, 

n = -1/g (||- + 1/2 v <f, . v <>) (17) 

By neglecting the last term, which is of the second order, 
the linearized equation for the free surface elevation is 

n = ~l/g §£ (18) 

On the surface, Z = 0, the linearized dynamic boun- 
dary condition (18) can be differentiated with respect to 
time and added to the linearized kinematic boundary con- 
dition (15) , to form a single boundary condition for the 
velocity potential: 



For plane progressive waves with finite-depth effects, 
the free surface elevation will be of the general form, 

n(x,t) = a cos (kx - wt + e) (20) 

in which the positive x - axis is in the direction of wave 
propagation, a is the wave amplitude, and e is some arbi- 
trary phase angle that can set equal to zero by the suit- 
able choice of the origin x = 0. The parameter k is the 
wave number, in which, 

k = 2 tt/ l 

where 

L = wavelength 

The solution to this problem is written in terms of 
a two dimensional velocity potential <|>(x,z,t), which must 
satisfy, the free-surface condition (19), Laplace's equa- 
tion, the bottom boundary condition S /S z = 0 at z = -d 
(for finite depth) , and must produce the wave elevation 
(20) for equation (18) . The potential must be sinusoidal 
due to (20) and the solution should be in the form: 

<f>(x,z,t) = Re[Z(z)e -ikx + iu3t ] (21) 

From Laplace's equation, Z must satisfy the ordinary dif- 
ferential equation: 



The most general solution to equation (22) 
of exponential functions is 

Z = Ce kz + De" kz 



in terms 



(23) 



where 

C = constant 
D = constant 

First assume the fluid depth is infinite or 
< z < 0. To avoid an unbounded motion deep below the 
free surface, the constant D in equation (23) must vanish 
leaving 

Z = Ce kz ' (24) 



If equation (24) is substituted into equation (21) , the 
velocity potential becomes, 



= Re[Ce k2 ’ ikx+i “ c ] 



(25) 



By setting z = 0, using equation (18) to calculate the 
linearized wave elevation, and comparing the result with 
the required plane wave in equation (20) , the constant 
C has the value of 



C = iga/u (26) 

For a fluid having a constant depth d, with the bot- 
tom being an impermeable plane, as stated earlier, the 
bottom boundary condition becomes, 



9 <J> 

9"z 



0 at z = -d 



(27) 



By looking at the general solution (equation (23)), it is 
found that both exponential functions must be retained, 
with constants C and D chosen such that equation (27) will 
be satisfied. 

For simplification, the two independent solutions 
e — are replaced by linear combinations of the hyperbolic 
sine and cosine function. For a fluid having finite depth, 
the solution becomes, 



Z(z) 



r cosh k (z + d) 
cosh kd 



(28) 



The constant C is selected such that equation (26) remains 
valid. When substituting (28) into (21) and using equa- 
tion (26) the velocity potential for a finite depth fluid 



becomes , 

ga cosh k(z + d) sin (kx _ ut) (29) 

Y 2 a) cosh kd 

Substituting equation (29) and (13) into equation 
(12), the combined wave-current velocity potential becomes: 

♦ - v + C °cosh ( kd ~~ ) - •"=> (30) 



Differentiating equation (30), the result becomes: 

U = = n + cosh k i^ Z i j- d ^- cos. (kx - cot) . . . . (31) 

9x c to cosh kd 



where 



U = combined wave-current velocity 

It can, therefore, be concluded that the combined 
wave-current velocity (U) can be obtained by simply add- 
ing the value for the uniform current velocity to the 
wave velocity in still water found using the parameters 
stated earlier. 

It is appropriate at this time to mention several 
points. Firstly, this relationship is invalid for non- 
linear boundary conditions. Also, since the wave pro- 
pagates into the flow, there will be modifications to 
both the waves and currents by interaction, therefore, 
the simple addition of current velocity and wave velocity 
is only approximate. Finally, linear wave theory is only 
valid for small amplitude waves, and many cases of wave- 
current interaction exist beyond the accepted range of 
validity for this theory. 

Wang (26) calculated the combined wave-current velo- 
city using the procedure outlined above and compared the 
calculated values with measured values. He derived a 
combined wave-current velocity relation given by: 

U = U c + ttH (1/T - U c /L) cosh [TT < z + d > 1 < 32 > 

1 ” 

where 

U = combined wave-current velocity 

U = uniform current velocity 

c 

By assuming that the distribution of current is uni- 
form with depth, the depth-averaged velocity was obtained 



from the integral of equation (32) . When comparing the 
calculated and measured combined wave-current velocities, 
a correlation coefficient of 0.939 was obtained, thus in- 
dicating this method for calculating currents is reasonabl 

Dimensional Analysis 

There are many variables encountered when studying 
scour. The following list is believed to be the important 
variables influencing scour depth at a cylindrical pile 
due to current and random wave motion. 



Description 


Symbol 


Units 


Dimensions 


Water depth 


h 


ft 


L 


Significant wave 
height 


H l/3 


ft 


L 


Kinematic viscosity 
of water 


V 


ft 2 /s 


l 2 /t 


Pile diameter 


D 


ft 


L 


Mean diameter 
of sediment 


D 50 


ft 


L 


Density of water 


P 


slugs/ft 2 


M/L 3 


Sediment/ fluid 
density difference 


< p s ■ p) 


slugs/ ft 2 


M/L 3 


Acceleration due 
to gravity 


g 


ft/ sec 2 


L/T 2 


Combined wave- 
current velocity 


U 


ft/ sec 


L/T 


Ultimate scour 
depth 


S d 


ft 


L 



Since random waves were used in this study, and they 
are defined by the significant wave height (H-^) , wave 



00 



period and wave length have been omitted from the above 
list of variables. 

The dimensional analysis was performed using a me- 
thod outlined by Venezian (25) during the 1984 Texas A&M 

University dredging short course. was chosen as the 

2 3 

unit of length, S d /v as the unit of time, and p/S d as 

the unit of mass. The remaining seven variables expressed 

in these units give the following dimensional numbers: 

S^/h , 3 , S^/D , S^D^q, p/ (p $ “ p) » v /^ d §> V /S^u 

By rearranging and combining terms : 

1 = UD = 

(S d /D) (v/S d u) v ln RP 

[(S d /h)(l/(v/S d u) 2 )(v 2 /S d 3 g)] % = u/^gh = N f 
[(S d /D)(l/(v/S d u) 2 )(v 2 /S d 3 g)] % = u/YiD = N Fp 
t(%)(v 2 /S d 3 g)(l/(v/S d u) 2 )] % = u/V2S^i= N es 

[(l/(v/S d u) 2 )(v 2 /S d 3 g)(S d /D 50 )(p/(p s - P))] % 

= U/V S D 50^ p s " p ) /p ) = N S 

Therefore, the functional relationship for scour is: 



-ff = f ( S d /H l/3’ S d /D ’ S d /D 50’ N RP’ N F’ N FP’ N ES ’ N S ) ' 



PART IV 



EXPERIMENTAL EQUIPMENT 

The scour experiments were conducted in a 120-foot 
long, 3-foot deep and 2-foot wide two dimensional, glass- 
walled wave flume. The test area was located approximate- 
ly 40 feet from the wave generator. The entire test area 
measured 40 feet in length and was divided into several 
sections. An 8-foot long 5-inch deep test sediment bed 
was bordered on each end by an 8-foot long 10- inch high 
false bottom. At the leading edge of the test area, a 
16-foot long ramp was used to gradually bring the wave 
up to the new depth at the top of the false bottom. A 
schematic of the test area is shown in Figure (22) . A 
1.5-inch aluminum pile was positioned along the center- 
line of the sediment bed. 

Currents were produced using a centrifugal pumping 
system. The flow rate was monitored using a Fischer and 
Porter electromagnetic flow meter. The pump discharge 
entered the wave flume just forward of the wave generator, 
and was diffused using a louvered grate. A Thermo-Systems 
model 1050 hot film anemometer was used to measure the 
vertical velocity distribution of the combined waves 
and currents . 

The random wave generator used was a Seasim modular 
wave making system. This system consisted of a servo 
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control system amplifier, a Seasim programmable spectrum 
random signal generator, a Seasim autocompensating wave 
height gauge, and a paddle type sea wave simulation roll- 
ing seal modular wave maker. A capacitance- type wave gauge 
was used with a Hewlett-Packard model 17403A carrier pre- 
amplifier to measure the wave profiles. To determine the 
significant wave height from the wave profile, a wave 
digitizer was utilized. 

A wave absorber was positioned in the wave flume at 
the opposite end of the wave generator to reduce any re- 
flections which might occur. 

Scour depth readings were made using a point gauge 



after each test. 



PART V 



EXPERIMENTAL PROCEDURES 

A total of 20 experiments were performed during this 
investigation. Four runs utilized currents only, and the 
remaining 16 runs used a combination of waves and cur- 
rents. All of the test runs used a 1.5-inch diameter pile 
placed in a bed of glass microbeads. The glass microbeads 
had a mean diameter (D^q) of 0.1 mm and a specific gravity 
of 2.45. For both the current only tests and the combined 
wave-current tests two water depths and two different cur- 
rents at each depth were used. Also used were four dif- 
ferent random wave spectra for each depth/current combina- 
tion. The data on the four random wave spectra are con- 
tained in Appendix A, and Table 1 contains a description 
of all tests conducted. 

At the commencement of each test, the sediment bed 
was leveled to a 5-inch uniform thickness. Water was 
added slowly to the tank so as not to disturb the sedi- 
ment. After the required water depth was attained, the 
proper flow rate through the tank was established. The 
water level in the tank was maintained by regulating a 10 
inch discharge valve located at the rear of the tank. 

Once the proper flow rate was attained and the required 
water depth held constant, the sediment bed was once again 
leveled. The signal generator, which had been previously 
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TABLE 1 - SUMMARY OF TESTS 



Test 



Current 



Wave Spectrum 



1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 



U 1 

U 2 

U 3 

U 4 

U 1 

U 2 

U 2 

U 2 

U 2 

U 1 

U 1 

U 1 

U 4 

U 4 

U 2 

U Z 

U, 

ul 

u~ 



D 

D 

P 

J1 

J2 

P 

J1 

J2 

D 

P 

J1 

J2 

D 

P 

J1 

J2 



U-^ = .405 fps 
U 2 = .380 fps 
= .650 fps 
= .485 fps 

D = Darbyshire 
P = Pierson - Moskowitz 
J1 = Jonswap (T = 8 sec.) 
J2 = Jonswap (T = 10 sec.) 



h 1 = 1.167 ft 
h 2 = .917 ft 







programmed, was started generating the desired wave spec- 
trum. The maximum scour depth was recorded at various in- 
tervals throughout the run, the time increment increasing 
as the test continued. Each test was conducted until an 
equilibrium in the measured maximum scour depth was at- 
tained. Following the completion of the test, the wave 
generator, and the flow through the tank were stopped and 
the water was slowly drained from the tank in order to 
preserve the scour pattern formed. Once the tank was com- 
pletely drained, photographs were taken of the scour pat- 
terns . In addition, sediment elevations, taken in the form 
of a rectangular grid, were measured using a point gauge 
in order to determine the equilibrium scour depth. This 
depth was used in the data calculations presented later 
in this paper. 



PART VI 



PRESENTATION AND DISCUSSION OF RESULTS 



The experimental data obtained from this study are 
presented below in Table 2: 

TABLE 2 - EXPERIMENTAL DATA 



Run 


Velocity (fps) 


Scour Depth (ft) 


H l/3 


1 


U c = 


.405 


.101 




2 


U = 
c 


.380 


.081 




3 


U C - 


.650 


.167 




4 


U C - 


.485 


.14 




5 


U = 


.442 


.126 


.2278 


6 


U = 


.394 


.115 


.1656 


7 


U = 


.393 


.13 


.1266 


8 


u = 


.401 


.127 


.1715 


9 


u = 


.386 


.09 


.1930 


10 


u = 


.437 


.122 


.1817 


11 


u = 


.459 


.152 


.1791 


12 


u = 


.419 


.113 


.2100 


13 


u = 


.498 


.143 


.1724 


14 


u = 


.491 


.149 


.1527 


15 


u = 


.504 


.174 


.1631 


16 


u = 


.508 


.155 


.1550 


17 


u = 


.665 


.18 


.1342 


18 


u = 


.659 


.174 


.0967 


19 


u = 


.660 


.176 


.1605 


20 


u = 


.664 


.168 


.1093 



The results of the laboratory experiments are present- 
ed in graphical form in this section. These graphs depict 
the relationship between the dimensionless parameters 



developed earlier for scour as well as several dimensional 
plots. On various dimensionless plots, data from experi- 
ments conducted by Wang (22) and Armbrust (2) were used 
in combination with data obtained from the random wave- 
current interaction experiment. On several plots, curves 
were constructed to indicate a trend only, and these plots 
were labeled accordingly. Regression analysis was per- 
formed on the remaining plots in order to construct the 
appropriate curves. Logarithmic regression analysis was 
used on the semi-logarithmic plots, and power regression 

2 

analysis was used on the full logarithmic plots. The "r " 
term adjacent to each curve represents the correlation 
coefficient for each curve. 

Figures (23) to (26) show the functional relationship 
that exists between the relative scour depth, S^/h, and 
S^/D and the sediment number for three distinct pile dia- 
meters and three different sediment sizes. It shows that 
scour depth increases with increasing wave-current velo- 
cities. Scour depth is not dependent on the sediment num- 
ber alone, but is also dependent on pile diameter (Figures 
23 and 25) and sediment size (Figures 24 and 26) . 

Since the basic mechanism for scour hole development 
is the vortex system resulting from flow separation, and 
the separation point is dependent on the Reynolds number, 
the strength of the vortex system is a function of the 
pile Reynolds number. Because of the relationship between 
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iig. 23. - Relative scour depth, S^/h, vs. sediment number for 
constant sediment size and various pile diameters 
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Fig. 24. - Relative scour depth, S^/h, vs. sediment number for 
constant pile diameter and various sediment sizes 
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26« — Relative scour depth, S^/D, vs. sediment number for 
constant pile diameter and various sediment sizes 



the vortex system and the pile Reynolds number, scour depth 
should, in turn, be functionally related to the pile Rey- 
nolds number. This functional relationship, along with 
Shen ' s curve for scour due to a steady current, is de- 
picted in Figure (27) . 

Figures (28) through (31) indicate the relationship 
between relative scour depth, S^/D and S^/h, and the pile 
Reynolds number. It shows that scour depth increases with 
increasing wave-current velocities or pile diameter as in 
the case of the sediment number; it can be seen from the 
graphs that scour depth is not dependent on the pile 
Reynolds number alone, but is also dependent on the pile 
diameter (Figures 28 and 30) and the sediment size (Figures 
29 and 31) . This agrees with the experimental work of 
Shen et al . (22); however, the dependence of the pile 

Reynolds number on sediment size differs from the conclu- 
sions reached by Palmer (18) in which he stated that the 
scour depth was independent of grain sizes from 0.12 mm 
to 0.63 mm. 

By comparing Figures (23) to (26) and Figures (28) to 
(31) with figures found in Lee (13) which plotted the same 
variables for multiple piles, it is noted that the single 
pile condition requires a higher velocity to produce the 
same scour depth as the multiple pile case, assuming all 
other variables (i.e., D^Q,h) are equal. 
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Fig. 27. - Equilibrium scour depth vs. pier Reynolds numb 
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Fig. 28. - Relative scour depth, S^/D, vs. pile Reynolds number for 
constant sediment size and various pile diameters 
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Fig. 29. - Relative scour depth, S^/D, vs. pile Reynolds number for 
constant pile diameter and various sediment sizes 
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Fig. 30. - Relative scour depth, S^/h, vs. pile Reynolds number for 
constant sediment size and various pile diameters 



lx 10' 



74 




m 

o 



X 



O 



x 




o 



b 



X 



X 







“O 

CO 



Fig. 31 • - Relative scour depth, S^/h, vs. pile Reynolds number for 
constant pile diameter and various sediment sizes 



It was found that if the pile Reynolds number and 
the sediment number were multiplied together and plotted 
against the relative scour depth, S^/h, all the points 
collapsed onto a single curve regardless of variance in 
pile diameter or mean sediment size. This result is 
shown in Figures (32) , (33) and (34) . The equation for 

each curve is also indicated on each figure. Due to the 
tendency of the data for random waves and currents and 
monochromatic waves and currents to collapse onto a sin- 
gle curve, it appears as if the irregularity of the waves 
is insignificant in terms of scour depth prediction. By 
using a dimensionless plot, as depicted in Figures (32)- 
(34) , it appears as though predictions of scour depth 
could be made for situations in which the pile diameter 
or mean sediment size varied without utilizing scaling 
factors. Before any conclusions can be made, field data 
are required to evaluate the need, if any, for scale fac- 
tors. 

Figures (35) through (37) show the relationship found 
for random wave-current interaction between the relative 
scour depth an< ^ the sediment number, N g , pile 

Reynolds number, Npp, and (N g ) (N^p) , respectively. The 
curves indicate that there is a rapid initial increase of 
the relative scour from the point of incipient motion up 
to a maximum relative scour depth. A further increase of 
N , N OT3 or (N ) (N dd ) will cause a rapid decrease in the 

S Kr S Kir 




(N, P )(N S ) 

Fig. 32. - Relative scour depth, S^/h, vs. sediment 

number times pile Reynolds number for a constant 
pile diameter and various sediment sizes 
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Fig. 33 . - Relative scour depth, S d /h, vs. sediment number times pile 
Reynolds number for a constant sediment size and various 
pile diameters 
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Fig. 34. - Relative scour depth, S d /h, vs. sediment' number times pile 

Reynolds number for various sediment sizes and pile diameters 
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Fig. 35 . - Relative scour depth, S d /H 1/3> vs< se diment number 
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Fig. 36. - Relative scour depth, S^/H, , , vs. pile Reynolds number 
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Fig. 37. - Relative scour depth, S./H-, ,? y vs. pile Reynolds number times 



relative scour depth, this decrease being independent of 
N s . or (N g ) (N^p) . These results are similar to those 

found by Wells and Sorenson (27) and presented earlier in 
this report in Figures (16) and (17) . The relative scour 
depth is smaller in Wells and Sorenson's case, but this 
can be explained by the fact that they used wave action 
only and, therefore, their scour velocities were smaller 
in magnitude. The pile Reynolds number for the rapid de- 

o 

crease of relative scour depth is approximately 7.5 x 10 
in both cases, this rapid decrease in S^/H-^ being indepen- 
dent of the pile Reynolds number and the sediment size. 

A dimensional plot was constructed in Figure (38) 
showing the relationship of ultimate scour depth, S^, ver- 
sus the combined wave-current velocity, U. It was found 
that in all cases the equilibrium scour hole was deeper 
than that associated with a steady boundary layer current 
with the same mean velocity. These results agree with 
predictions made by Niedoroda and Dalton (15) . In gener- 
al, it can be stated that -the scour depth caused by wave- 
current interaction is approximately 10% greater than the 
scour depth due to a steady current alone. 

In Figure (39), scour depth was compared to H-^^. 
Although the data are limited, it appears as though for 
current velocities greater than 0.45 feet per second, there 
is a tendency for scour depth to increase with increasing 
wave height, for velocities less than 0.45 feet per second, 
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Fig. 38. - Ultimate scour depth vs. current velocity and combined 
current-wave velocity 
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Fig* 39* - Ultimate scour depth vs, significant wave height 



scour depth decreases with increasing wave height and for 
velocities approximately equal to 0.45 feet per second, 
there is no definite trend at all. This result is most 
likely due to the fact that waves and currents interact 
differently as the magnitude of the current changes. Al- 
though a discussion was presented earlier in this paper 
on the mechanics of scour due to waves and currents, there 
is still a great deal of uncertainty in this area, and 
further research is indicated. 

Scour Pattern Observations 

The rate of scour hole development by wave-current 
interaction was faster than that due to a steady current 
alone. A typical example of scour development is shown 
in Figure (40) . As seen in this figure, scour development 
is very rapid initially, the scour rate then decreases 
until the equilibrium scour depth is reached. 

In all cases, the scour hole due to wave-current in- 
teraction was similar in shape to the scour hole formed by 
a steady current. This shape resembles an inverted cone. 

A typical example of a scour hole formed by a steady cur- 
rent is shown in Figure (41) . Figures (42) to (45) depict 
the scour pattern formed by the same current used in Run 
4 (Figure 41) , combined with four different random wave 
spectra. Photographs of the scour patterns formed by 
wave-current interactions for the remaining 12 runs can 
be found in Appendix B. As stated earlier, the equilibrium 
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Fig. 40. - Ultimate scour depth vs. time 



87 




Fig. 41 . - Typical scour cavity due to a 
steady current 



h = 1.167 ft 
U = .485 ft/sec 
D c = .125 ft 
S d = .14 ft 
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Fig. 42 . - Scour pattern due to random waves and 
current - Run 13 . 



h = 1.167 ft 
U = .485 ft/sec 
IT = .498 ft/sec 
wave spectrum = Darbyshire 

h 1/3 = .1724 ft 

D = .125 ft 

s d = .143 ft 




Fig. 43. - Scour pattern due to random waves and current 
Run 14. 



h = 1.167 ft 
U = .485 ft/sec 
U c = .491 ft/sec 
wave spectrum = Pierson' 

H 1/3 = .1527 ft 

D = .125 ft 
S d = .149 ft 



Moskowitz 
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Fig. 44 . - Scour pattern due to random waves 
and current - Run 15. 



h = 1 .167 ft 
U = .435 ft/sec 
U c = .504 ft/sec 
wave spectrum = Jonswap (T 

H 1/3 = .1631 ft 

D = .125 ft 

s d = .174 ft 



= 8 sec 
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Fig. 45. 



Scour pattern due to random waves 
and current - Run 16. 



h = 1.167 ft 
U c = .485 ft 

U = .508 ft 

wave spectrum = Jonswap (T = 10 se 
H 1/3 = .1550 ft 

D = .125 ft 
S d = .155 ft 
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scour hole associated with the wave-current interaction 
was deeper when compared to the equilibrium scour hole 
formed by a steady current alone. The size and shape of 
the equilibrium scour hole associated with wave-current 
interaction was at times greater than, approximately equal 
to, or less than the size and shape of the scour hole due 
to a steady current alone. These results differ from the 
prediction made by Niedoroda and Dalton (15) . They stated 
that the scour hole size and shape for combined waves and 
currents should not be as great as that developed by a 
steady current alone. 

Scour development is illustrated by Figures (46) to 
(48) . Figure (46) shows a typical pattern of scour develop- 
ment in its initial stages, approximately the first 2 to 
5 minutes of the test. Figure (47) and (48) are the typi- 
cal scour patterns that occur due to wave-current interac- 
tion 30 minutes and 2 to 3 hours after test commencement, 
respectively . 

The general scour pattern was one of localized scour 
around the pile. A vortex formed around the pile due to 
disturbed flow and suspended the sediment. The sediment 
was carried farther downstream under the combined effects 
of waves and currents than it was with a steady current 
alone. Initially, the sediment suspension appeared to 
become greater as the scour hole developed, possibily due 
to a strengthening of the vortex system surrounding the 
pile . 
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Fig. 46 



Initial stages of scour development 





yn- 




Fig. A- 7. 



Scour pattern 30 minutes after 
test commencement 
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Fig. 48. 



Scour pattern 2 to 
test commencement 



3 hours after 
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Because of the disturbance of flow caused by the pile, 
trailing ripples formed immediately behind the pile ex- 
tending in the direction of wave propagation as the test 
progressed. In nearly all cases, fairly uniform ripples 
formed along the entire sediment bed. The wave length of 
these ripples ranged from approximately 1.5 to 2 inches. 

Most of the ripples that covered the sediment bed initial- 
ly formed at the interface between the false bottom upstream 
of the pile, and the sediment bed. The ripples probably 
formed at this location due to a disturbance in flow, and 
once formed, the ripples progressed in the direction of 
wave propagation. 

In all cases, the maximum scour depth occurred down- 
stream of the pile, approximately 1/4 inch from the pile, 
and at an angle of approximately 40° to the left or right 
of an imaginary line drawn through the pile and parallel 
with the wave flume walls. 
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PART VII 
CONCLUSIONS 

1. Scour development was very rapid initially. The scour 
rate then decreased until an equilibrium scour depth 
was reached. 

2. Scour development due to current and wave action oc- 
curred faster than scour due to a steady current alone. 

3 . The equilibrium scour depth caused by wave-current 
interaction is approximately 107 o greater than the scour 
depth caused by a steady current alone. 

4. The shape of the scour hole due to waves and currents 
and the shape of the hole due to current alone are 
similar. The shape resembles an inverted cone. 

5. For a constant pile diameter and sediment size, scour 
depth increases with increasing wave-current velocity. 

6. For a constant wave-current velocity and sediment size, 
scour depth increases with an increase in pile diameter. 

7. For a current velocity greater than 0.45 feet/ sec, 
scour depth increased with increasing wave height. 

For a current velocity less than 0.45 feet/ sec, scour 
depth decreased with increasing wave height. 

8. A relationship exists between scour depth and sediment 

number (N ) . However, scour depth is not dependent on 
s 

N g alone, but is also dependent on pile diameter and 



sediment size. 



9 . 



A relationship also exists between scour depth and 
the pile Reynolds number (N^p) . Again, scour depth 
is not dependent on N^p alone, but is also dependent 
on pile diameter, as well as sediment size. 

10. By plotting the relative scour depth (S^/h) against 
the product of the pile Reynolds number and sediment 
number , (Npp) (N g ) , all data points collapse onto a 
single curve regardless of variance in pile diameter 
or sediment size. 

11. A relationship exists between the relative scour depth, 

(S d /H l/3>’ and N RP’ N s and (N s )(N RP } ' 

12. In order to attain the same scour depth with a single 

pile as that achieved in a multiple pile case, the 
velocity must be greater assuming all other variables 
are identical. 

* 1 

Although the amount of data is limited, due to the 
fact that all data could be collapsed onto one curve when 
plotting S^/h versus (N^p) (N g ) , (Figures 32 and 34) , it 
appears as though the type of wave train used for model 
studies, either monochromatic or random, is insignificant 
in the prediction of scour depth about a single cylindri- 
cal pile. However, when attempting to predict the geometry 
of the scour hole, the irregularity of the wave train be- 
comes very significant. Lee (13) and Armbrust (2) stated 
that the pattern of scour for wave-current interaction 



is no longer conical in shape, but the pattern is irregular 
due to the wave particle orbital motion. The data in this 
study clearly indicate a conical scour pattern for random 
wave-current interaction. Although these results are dif- 
ferent from those of Lee and Armbrust, they agree with 
field data, as well as with the predictions of Niedoroda 
and Dalton (15) . The fact that the geometric shape of 
the scour hole differs with the type of wave train used, 
may become important in the prediction of scour depth 
around multiple pile cases . 

Figures (35) and (36) show the relationship found for 
random wave-current interaction between the relative scour 
depth an <3 the sediment number, N g , and pile Rey- 

nolds number, N^p, respectively. These results are simi- 
lar to those found by Wells and Sorenson (27) and present- 
ed in Figures (16) and (17) . Because of similar trends 
in all figures and due to the fact that the pile Reynolds 

number for the rapid decrease of relative scour depth is 

3 

approximately 7.5 x 10 in both cases, (Figures 36 and 17), 
it again appears as though the irregularity of the wave 
train is insignificant . in the prediction of scour depth 
around a single cylindrical pile. 

In addition, Armbrust (2) found that there was an in- 
crease of scour depth with increasing wave height. As 
stated earlier, for irregular waves when the current ve- 
locity was greater than .45 ft/ sec, scour depth increased 



with increasing wave height, and the scour depth decreased 
with increasing wave height for a current velocity less 
than .45 ft/sec. 

The conclusions resulting from this study were found 
for a flow in which the steady current was greater than 
the oscillatory flow produced by the waves. If the mag- 
nitude of the two flows in relation to one another change, 
the flow patterns will also change and the geometry of 
the scour hole may very likely differ from results found 
in this study. 

As mentioned in previous investigations, to attempt 
to predict scour depth for a prototype using the results 
of this study could lead to erroneous results. In order 
to accurately predict the prototype behavior, complete 
similitude is required between model and prototype, and 
field data is required to evaluate scale effects, should 
any exist. 



Suggestions for Further Research 

1. This study could be repeated under the same conditions, 
thereby omitting scale factors, if any exist, but using 
monochromatic waves to determine if the irregularity 

of the wave train is significant in scour studies. 

2. Further research can be conducted to investigate the 

variation of scour depth with varying significant wave 
height, as the current velocity varies. 
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3. It would be desirable to investigate why there is a 

rapid decrease in relative scour, f° r a pile 

3 

Reynolds number, N^p, of approximately 7.5 x 10 . 

4. The study could be repeated using various sediment 
sizes and pile diameters to determine if scale effects 
are negligible when using a plot of S^/h versus (N^p) 
(N s ). 

5. This study could be repeated using multiple piles in 
order to compare the results with previous studies in- 
vestigating scour about multiple piles due to monochro- 
matic wave-current interaction. 
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APPENDIX A 



Random Wave Spectra Data 
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iTER DEPTH = .61 m 

JQNSWAP SPECTRUM 

3 = 10.4 , T = 10 , F0 = . 1 , Y = 3.3 , calculated HS = 10. 267592SS55 



[SPLACEMENT CONTROL 



WE WAVE TRANSFER FUNCTIONS SIGNAL FILTER MODEL 



iEQUENCY 


SPECTRUM 


TRANS 1 


TRANS2 


FINAL 


SPECTRUM 


FUNCTION 


FREQUENC't 


Hz 


m'2/Hz 








mm--2/Hz 




Hz 


.01313 


0.0000 


9. 6217 


1 . 0343 


9.9513 


0.00 


0 . 0000 


. 131 : 


. 02625 


0 . 0000 


4. 6596 


1 . 0635 


4.9738 


0 . 00 


0.0000 


. 262= 


. 039w3 


0 . 0000 


2. 9445 


1. 1028 


3. 2471 


0. 00 


0. 0000 


.3936 


. 05250 


. 0004 


2.0507 


1. 1370 


2. 3317 


.02 


. 0068 


. 525C 


. Uo56o 


nt nng^i 


1 . 5002 


1. 1713 


1 . 7572 


63 . 80 


. 3794 


. 656" 


.07875 


29.6923 


1 . 1417 


1 . 2055 


1 . 3764 


562. 52 


1 . 0397 


.737= 


.09188 


111.3599 


' .9137 


1 . 2398 


1. 1327 


1428.88 


1 . 7460 


.9186 


. 10500 


17S. 6736 


.7776 


1 . 274 1 


. 9907 


1753.74 


1.9461 


1 . 050C 


. 11313 


61.3725 


. 6991 


1 . j08o 


.9146 


513. 40 


1 . 0602 


1. 18 IT; 


. 13125 


38.7328 


. 6516 


1.3426 


. 8747 


296.38 


.8118 


1 . 3125 


. 14438 


27. 4067 


. 6201 


1 . 3768 


. 8538 


199. 77 


.6723 


1 . 4436 


. 15750 


19.3036 


.5978 


1.4111 


. 8435 


137.34 


.5627 


1 . 575C 


. 17063 


13. 6772 


.5812 


1 . 4453 


. 8400 


96.51 


.4766 


1 . 706z 


. 13375 


9 . 3063 


.5635 


1.4796 


.3412 


69.39 


. 4087 


1.8375 


. 19688 


7. 1314 


.5587 


1 . 5138 


. 8458 


51.01 


.3546 


1.9686 


. 21000 


5. 2632 


. 5509 


1 . 5481 


.3523 


38.28 


.3112 


2. 1 00 C 



rDEV PADDLE DISPLMT = 26.1168968025 mm 

*X STROKE = 129.659502619 mm 4.96458302834 STDEVS 

LX SIGNAL OUTPUT = 7.09289198229 STDEVS 
LIGHTING FUNCTION 



+ 0. 00000 
-6.97203 
-8.77976 
+6. 01734 
-4.84562 
+ . 77440 
- 1.01222 
+ . 36305 

:x VALUE STARTING AT ADDRESS OOCO 

> 72 7F 35 ED DO CC D1 E8 

> 00 07 05 02 01 FC F8 FD 



+41 . 85458 
-17.37474 
+4. 30239 
+2. 03808 
-. 15382 
+. 18659 
+ 1 . 93385 
-. 19043 



+46.5143 
-18.9913 
+11. 9233 
-2.5298 
+2. 3983 



OC 

05 



- 1 . 
+ 2 . 
-1 . 

21 

07 



5718 

4316 

7046 

ID 

03 



.OCK INTERVAL = .119047619043 S = ID HEX 
>T I MATED GAIN SETTING = 6.33577130898 



10 

01 



+19.5321 
-17. 2952 
+10. 5649 
-6. 1335 
+1.6312 
-2.8261 
+1.0725 
-2. 1668 

06 F9 EF 
FF FB FA 



OB 02 OF 01 11 03 

2K 32K 1 28K 



lEDBACK OPTIONS: 
: eedback Points OA 03 
sequence Length IK 



14 03 
1M 



vG CD xl LI t-J 1- O -0 '-J O LI L4 10 

O' it O *4 *-* CD LH h- 00 LI 10 -O O' it 

m 'O O' tn it t o h o co o a- li tt t o 

m f n t i i m r n i .1 m in f 1 r- < m r n t .1 
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'* 5>int-ui_3E. ruNL, i xuin *** 
iDEL SCALE = Is 100 
ITER DEPTH = . 61 m 

JONSWAP SPECTRUM 

| = 10.4 , T » 8 , FO = .125 , Y = 3.3 , calculated HS = 10.1079609998 



3PLACEMENT CONTROL 

7E WAVE TRANSFER FUNCTIONS SIGNAL FILTER MODEL 

IQUENCY SPECTRUM TRANS 1 TRANS2 FINAL SPECTRUM FUNCTION FREQUENCY 



Hz 


m"'2/Hz 






mm'‘-2/Hz 






Hz 


. 0 1 j 1 o 


0 . 0000 


9.6217 1.0343 


9.9513 


0. 00 


0 . 0000 




. 1313 


. 02625 


0.0000 


4.6596 1.0685 


4.9788 


0.00 


0. 0000 




.2625 


. 03938 


0. 0000 


2.9445 1.1023 


3.2471 


0 . 00 


0 . 0000 




.3938 


. 05250 


0.0000 


2.0507 1.1370 


2. 3317 


0 . 00 


0.0000 




. 5250 


. 06563 


. 0003 


1 . 5002 1.1713 


1 . 7572 


. 01 


. 0046 




. 6563 


. 07875 


.6722 


1.1417 1.2055 


1 . 3764 


12.73 


. 1640 




.7875 


.09188 


1 1 . 9993 


.9137 1.2398 


1 . 1327 


153.97 


.5731 




.9188 


. 10500 


39.5000 


.7776 1.2741 


. 9907 


337.70 


.9150 


1 


. 0500 


. 11813 


123. 9536 


.6991 1.3083 


.9146 


1036.90 


1 . 5068 


1 


. 1813 


. 13125 


145.2539 


.6516 1.3426 


.8747 


1111.46 


1.5721 


1 


.3125 


. 1 4438 


58.8910 


.6201 1.3768 


. 8538 


429.27 


.9355 


1 


.4438 


. 15750 


36 . 3984 


.5978 1.4111 


.8435 


258.96 


.7727 


1 


. 5750 


. 17063 


27.4429 


.5812 1.4453 


. 8400 


193. 63 


. 6751 


1 


. 7063 


. 18375 


20.7723 


. 5685 1 . 4796 


.8412 


146.99 


.5948 


1 


. 8375 


. 19688 


15. 6930 


.5587 1.5138 


.8458 


112.26 


.5261 


1 


. 9688 


.21000 


11.9024 


.5509 1.5431 


. 8528 


86.56 


. 4680 


*-y 


. 1000 


DEV PADDLE DISPLMT = 


22.5873917467 mm 










< STROKE 


= 129.659502619 mm 5.74034860122 STDEVS 










< SIGNAL 


OUTPUT = 7. 


38692477717 STDEVS 










IGHTING 


FUNCTION 
















+0.00000 


+41.74525 


+37.5259 


+. 2737 










-24. 19653 


-19.74887 


-6.0267 


+2. 0089 










+7. 45484 


+10.32406 


+5.6980 


-4.5183 










-7.68504 


-2.01908 


+2.6540 


+1.9774 










+.57594 


+.24903 


- 1 . 5343 


-3.0336 










-.40299 


+3. 45594 


+2.8269 


-. 7413 










-1.74790 


-.47922 


-. 5562 


-1.0337 










+ . 56383 


+1 . 89930 


+.2708 


-1 . 3650 








( VALUE 


STARTING AT 


ADDRESS OOCO 












7F 72 


01 B6 C4 


EE 06 17 21 


11 F2 E9 


FA 08 06 








01 FB 


F7 FF OB 


09 FE FB FF 


FE FD 02 


06 01 FC 









DCK INTERVAL = .119047619048 S = ID HEX 
r I MATED GAIN SETTING = 5.70669069914 



ID BACK OPTIONS: 

sedback Points OA 03. OB 02 OF 01 11 03 14 03 

rqusnce Length IK 2K 32K 128K 1M 



-OCITY CONTROL 
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t* S IMPULSE FUNCTION *** 

1DEL SCALE - 1: 100 

4 TER DEPTH = .61 m 

PIERSON SPECTRUM 

4 = 10 , T = 10 , FO = .0719942404608 . Y = 1 , calculated HS = 

.91768977232 



JSPLACEMENT CONTROL 

I YE WAVE TRANSFER FUNCTIONS SIGNAL FILTER MODEL 



'EQUENCY 


SPECTRUM TRANS 1 


TRANS2 


FINAL SPECTRUM FUNCTION 


FREQUENCY 


He 


(tK-2/H= 








mm "-2 /Hz 




Hz 


.01313 


0. 0000 


9. 6217 


1 . 0343 


9. 9513 


0. 00 


0 . 0000 


■ 1 1 O 


. 02625 


0. 0000 


4. 6596 


1 . 0685 


4.9788 


0.00 


0 . 0000 


. 2625 


. 03938 


.0149 


2.9445 


1. 1028 


3. 2471 


1. 57 


. 0570 


.3933 


. 05250 


30.2292 


2.0507 


1. 1370 


2.3317 


1643. 48 


1 . 8477 


. 5250 


. 06563 


117.5497 


1 . 5002 


1. 1713 


1 . 7572 


3629. 65 


2.7559 


. 6563 


. 07875 


114.8723 


1. 1417 


1.2055 


1 . 3764 


2176.26 


2. 1434 


.7875 


.09188 


77.7870 


.9137 


1 . 2398 


1. 1327 . 


998. 10 


1 . 4592 


.9188 


. 10500 


47.9994 


. 7776 


1.2741 


. 9907 


471 . 13 


1 . 0087 


1 . 0500 


. 1 1813 


29.3902 


.6991 


1 . 3083 


. 9146 


245.86 


. 7337 


1.1813 


. 13125 


18. 3583 


. 65 1 6 


1 . 3426 


. 8747 


140. 47 


.5589 


1.31 25 


. 14438 


11.7933 


. 6201 


1 . 3768 


. 8538 


85.96 


.4410 


1 . 4438 


. 15750 


7.7991 


. 5978 


1.4111 


• 84o5 


55.49 


.3577 


1 . 5750 


. 17063 


5.3014 


. 5812 


1.4453 


. 8400 


37.41 


.2967 


1 . 7063 


• 1 So / 5 


3.6953 


. 5685 


1 . 4796 


.8412. 


26. 15 


. 2509 


1 . 8375 


.19688 


2.6349 


.5587 


1.5138 


.8458 


18. 85 


.2156 


1.9688 


.21000 


1 . 9174 


. 5509 


1 . 5431 


.8523 


13.94 


. 1378 


2. 1000 


DEV PADDLE DISPLMT = 


35.3804787138 mm 








X STROKE 


= 129.659502619 mm 3.66471871869 STDEVS 








X SIGNAL 


OUTPUT = 6. 


51233980254 STDEVS 








RNING, TRUNC WAVES 














IGHTING 


FUNCTION 
















+0.00000 


+42. 


63063 


+57.7463 


+45.7075 








+25. 71290 


+9. 


1 1734 


-5. 5703 


-18. 4574 








-24.36528 


^ 


24950 


-17. 0950 


-11.7322 








-4.92107 


+2. 


56339 


+7.2982 


+8. 6747 








+9. 12731 


+9. 


0933 1 


+7. 0203 


+3.6188 








+1 . 26716 


-. 


0 1 803 


-1 . 7862 


-3 . 50 1 9 








-3. 60015 




79689 


-2.7233 


-2. 8087 








-1.81417 


^ • 


65941 


-.6409 


-.8160 






X VALUE 


STARTING AT 


ADDRESS 


OOCO 










5E 7F 


65 39 1 4 


F4 D7 CA CF 


DA E6 F5 


06 10 13 






14 OF 


08 03 00 


FC FS F8 FA 


FA FA FC 


FF FF FE 






QCK INTERVAL = .119047619048 S = ID HEX 









TIMATED GAIN SETTING = 7.38052823039 
EDBACK OPTIONS: 

esdback Paints OA 03 OB 02 OF 01 11 03 14 03 

equence Length IK 2K 32K 128K 1M 



LOCI TV CONTROL 



t* S IMPULSE FUNCTION *** log 

IDDEL SCALE = 1: 100 

i^TER DEPTH = .61 m 

DARBYSHIRE SPECTRUM 

Is = 10.4 , T = 10 , FO = . 1 Y = 1 , calculated HS = 9.98580007236 



^PLACEMENT CONTROL 

UVE WAVE TRANSFER FUNCTIONS SIGNAL FILTER MODEL 



JEOUENCY 


SPECTRUM 


TRANS 1 


TRANS2 


FINAL 


SPECTRUM 


FUNCTION 


FREQUENCY 


Hz 


m-"-2/Hz 








mm-'-2/Hz 




Hz 


.01313 


0 . 0000 


9. 6217 


1 . 0343 


9.9513 


0 . 00 


0.0000 


. 1313 


. 02625 


0.0000 


4.6596 


1 . 0685 


4.9788 


0 . 00 


0 . 0000 


.2625 


. 03938 


0. 0000 


2.9445 


1 . 1028 


3. 2471 


0 . 00 


0 . 0000 


. 3938 


. 05250 


0 . 0000 


2. 0507 


1 . 1370 


2.3317 


0 . 00 


0 . 0000 


. 5250 


. 06563 


2.2426 


1 . 5002 


1. 1713 


1 . 7572 


69.25 


. 3807 


. 6563 


. 07875 


32. 5300 


1. 1417 


1 . 2055 


1.3764 


616.28 


1. 1406 


.7875 


.09188 


100.0229 


. 9137 


1.2398 


1 . 1327 


1283. 41 


1.6547 


.9188 


. 10500 


125.7719 


.7776 


1.2741 


. 9907 


1234. 49 


1 . 6328 


1 . 0500 


. 1 1813 


72.3776 


. 6991 


1 . 3083 


.9146 


605.45 


1. 1514 


1 . 1313 


.13125 


46.0842 


. 6516 


1.3426 


. 8747 


352.63 


.8855 


1.3125 


. 14438 


31.3266 


. 620 1 


1.3768 


. 8538 


228. 35 


.7187 


1 . 4438 


. 15750 


22.2956 


. 5978 


1.4111 


. 8435 


158.62 


. 6047 


1 . 5750 


. 17063 


16.4187 


.5812 


1 . 4453 


. 8400 


115.85 


.5221 


1 . 7063 


. 1 3375 


12.4146 


. 5685 


1 . 4796 


.8412 


87.85 


.4598 


1 . 8375 


. 19688 


9.5875 


.5587 


1 . 5138 


.8458 


68. 58 


. 4112 


1.9688 


. 21000 


7 • 5 ■ j* c* 6 


. 5509 


1.5431 


.8528 


54.79 


.3723 


2. 1000 



f*DEV PADDLE DISPLMT = 25.2254078414 mm 

'iX STROKE = 129.659502619 mm 5.14003592362 STDEVS 

Kx SIGNAL OUTPUT = 7.31375797316 STDEVS 

LIGHTING FUNCTION 





+0. 00000 


+42.31362 


+46 . 0740 




+ 17. 


2287 




-8.53612 


-16.32912 


-16.5965 




-15. 


6599 




-8.43808 


+3. 58199 


+10. 0639 




+8. 


1905 




+4.94569 


+3. 34469 


+. 1 129 




—3. 


9888 




-4. 13193 


-.99501 


+.2543 






7096 




-.25819 


+ 1 . 1411 1 


+.3952 




- 1 . 


1932 




-. 35002 


+1 . 47241 


+1 . 0045 




“ * • 


4428 




-.05661 


+.71136 


-.4055 




- 1 . 


3833 


•:x VALUE 


STARTING AT 


ADDRESS OOCO 










> 76 7F 


2F E9 D3 


D2 D5 E9 


OA 1C 17 


OE 


09 


00 1 


i FD 0 1 


FE FF 03 


01 FD FF 


04 03 FF 


00 


02 


FF 1 


.OCK INTERVAL = .119047619048 S = 1 


D HEX 








5TI MATED 


GAIN SETTING 


= 6.31005805647 








iEDBACK OPTIONS: 












"esdbac k 


Foints OA 03 


OB 02 OF 01 


11 03 14 


03 






Sequence 


Length IK 


2K 32K 


1 28K 


1M 







iLOCITY CONTROL 



APPENDIX B 



Scour Patterns for Random Wave-Current Interaction 



Ill 




Run 5* 



h = 1.167 ft 
U c = .405 fps 
U = .442 fps 
wave spectrum = 
H y/y. ,227% ft 
D = .125 ft 
S d = .126 ft 



Darbyshire 
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Run 6. 



h 


11 

• 

'vD 


ft 


U c 


II 

• 

Oo 

O 


fps 


U 


= .394 


fps 



wave spectrum = Darbyshire 
H 1// ,= .16% ft 
D = .125 ft 
S d = . 1 1 5 f t 
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Run 7. 



h = .91? ft 

U c = .360 fps 

U = .393 fps 

wave spectrum = Pierson— Moskowitz 
H'j ^ 266 ft 

D = .123 ft 

S d = .13 ft 
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Run 8. 



h = .917 fp 

U c = .380 fps 
U = .4-01 fps 
wave spectrum = Jonswap (T 

h 1/3 = .1713 ft 

D = .125 ft 

s d = .127 ft 



= 8 sec’ 
< 
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Run 9. 



h = .917 ft 

U c = .380 fps 

U = .386 fps. 

wave spectrum = Jonswap (T = 10 s< 

h i/ 3= .1930 ft 
D = .125 ft 

s d = .09 ft 
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Run 10. 



h = 1 .167 ft 
U c = .405 fps 
u = .437 fps 

wave spectrum = Pierson-Moskowitz 
h 1/ 3 = .1317 ft 
D = .125 ft 
S d = .122 ft 
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Run 1 1 . 



h = 1.167 ft 
U c = .405 fps 
U = .459 fps 

wave spectrum = Jonswap (T = 8 sec 

H 1/3= * 1791 ft 
D = .175 ft 

S d = .152 ft 




118 



I 




Run 12. 



h = 1 .167 ft 

U c = ,/f 05 fps 

U = .419 fps 

wave spectrum = Jonswap (T 

H .210 ft 

D = .125 ft 

s d = .113 ft 



10 S€ 
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Run 17. 



h = .917 ft 

U c = .650 fps 

U = .665 fps 

wave spectrum = Darby shire 

h 1/3 = .1342 ft 

D = .125 ft 

S d = .18 ft 




Run 1 8* 



h = .917 ft 
U c = .650 fps 
U = .659 fps 
wave spectrum = 
E = .0967 ft 
D = .125 ft 

s d = .174 ft 



pier so n— Moskowit z 
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Run 19* 



h = .917 ft 
U c = .650 fps 
U = .660 

wave spectrum = Jonswap (T 
.1605 ft 
D = .125 ft 
S d = .176 ft 






= 8 se< 
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Run 20. 



h = .917 ft 
U c = .650 fps 
U = . 664 fps 

wave spectrum = Jonswap(T = 1 0 se 
H 1/3 = .1093 ft 
D = .125 ft 
S d = .168 ft 
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APPENDIX C 
List of Notations 



Symbol 


Quantity 


Units Dimensions 


a 


wave amplitude 


ft 


L 


D 


pile diameter 


ft 


L 


D 50 

% 


mean diameter of sediment 
acceleration due to gravity 


ft 

ft/ sec 2 


L 

L/T 2 


h 


still water depth 


ft 


L 


H l/3 


significant wave height 


ft 


L 


H 


wave height 


ft 


L 


L 


wave length 


ft 


L 


n es 


scour Euler number 




, , 


n f 


flow Froude number 






n fp 


pile Froude number 






n rp 


pile Reynolds number 






N 

s 


sediment number 






s d 


ultimate scour depth 


ft 


L 


T 


wave period 


sec 


T 


U 


mean wave- current velocity 


ft/ sec 


L/T 


u 

c 


current velocity 


ft/ sec 


L/t 


U 


Ursell parameter 


. - 


— 


IT 

V 


kinematic viscosity of water 


ft 2 /s 


l 2 /t 


P 


fluid density 


slug/ft 2 


M/L 3 


P c " P 


sediment /fluid density 


slug/ ft 2 


M/L 3 


s 


difference 




CO 


angular velocity of 
oscillation 


rad/ sec 


1/T 


0 


angle of repose 




r\ 


T 


average shear stress 


lb/ft 2 


M/Lx 


n 


water surface elevation 


ft 


L 



Thesis 
Ell 3 
c. 1 


z*';83^G 

Eadie 

Scour about a single , 
cylindrical pile due to 
random waves and cur- 
rent . 



Thesis 

E113 

c.l 


. , i * n f ^ 

Eadie 

Scour about a single , 
cylindrical pile due to 
random waves and cur- 
rent . 



thesE 113 

Scour about a single, cylindrical pile d 
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